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Question 1 [10 Marks] 
(a) A vector is input in MATLAB via the command: 
a= [1 49 16 25]; 


Write down the expected output from the following subsequent commands: 


(i) b = sqrt(a) 
(ii) c = a(a<i15) 
(iii) d = find(a>=8) 
(iv) e = [a;-a] 

(v) f= a°2 


(b) Consider the following function, func.m: 


function [f] = func(x,y) 
f=x"2-y"2; 


end 


What would be the associated output from the MATLAB command z=func (3,2)? 


(c) Assuming that x and y are vectors of real numbers, what does the following MATLAB 


function do? 


function m = func2(x,y) 


n = min(length(x) ,length(y)); 


for ii = 1i:n 
if x(ii) > y(ii) 
m = mt1; 
end 
end 


end 
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Question 2 [10 Marks] 


(a) Write a MATLAB function to determine the average of all the elements of a matrix 


mat without using MATLAB’s built-in mean command. 


(b) Write a MATLAB script to evaluate y = sinz and z = cosz at 1000 equally spaced 
values of x over the interval —7 < x < 7, and then graph the functions as a red dashed 
line and a blue solid line, respectively. Add code to your script to label both the x- 
and y-axes of the graph, to add a title and a legend. 
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Question 3 |4 Marks] 
(a) Convert the following numbers from binary to decimal form, showing all working: 
(i) 110, (iii) 110.101, 


(ii) 0.101, (iv) 101.0011, 


(b) Convert the following numbers from decimal to binary, using the pen and paper algo- 


rithms demonstrated in lectures: 
(i) 51 (ii) 0.25 (iii) 51.25 


(c) Re-write the binary representation of 51.25 as a double precision floating point number 
(in binary) according to IEEE standard 754. 
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Question 4 [3 Marks] A standard form for the Taylor series of a suitable function f(x) 


about a point x7 = a is: 


"a "a (n) 
g(x) = f(a) + HO @ — a) ¢ Fe aye t+ Ee ay" + Basal, 
where p(y.) 
= Ln —a n+1 
Rn4i(x) (n+ 1)! (@ ) ? 


and jl, is some number between x and a. Use the above form of the Taylor series to 
determine the first three non-zero terms of the Taylor series expansion of cosx about 


x = 7/2 (you are not required to expand any brackets to simplify your Taylor series). Use 
COs £ 


x—m/2 


the terms of your series to then find an approximation of when z is close to 77/2 


(again do not expand any brackets). 


AMTH250 Exam Trimester 1 2018 


Question 5 |7 Marks] 
The following is an incomplete section of MATLAB script to apply the bisection algorithm 


to find a zero of the function contained in a file func3.m within the interval [2, x,]. 


clear 
xl = input(‘Enter left end point of interval containing solution: ’); 


xr = input(‘Enter right end point of interval containing solution: ’); 


nmax = 30; 
fl = func3(x1) ; 
for ii = 1:nmax 
xnew = ... %Incomplete line 1 


fnew = ... 4Incomplete line 2 


if flxfnew > 0 
xl = xnew; 


fl 


... Alncomplete line 3 
else 
xr = ... 4Incomplete line 4 


end 


if abs(xr-xl) < ... %Incomplete line 5 
break 
end 


end 


(a) Complete the indicated lines of code so that the zero is determined correct to 8 decimal 


places. 


(b) What is the maximum length of the starting bracketing interval such that the script will 
actually terminate with a solution which is guaranteed to have the required accuracy? 
(Note the line nmax = 30;) 
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Question 6 [10 Marks] 


Consider the function g(x) = cos x. 


(a) 


Sketch a simple graph of y = g(x) over the interval 0 < x < 27, along with a graph of 
y = x. Based on your graph, identify an approximate solution to x = g(x) that could 
be used as the initial input, xo, into a fixed point iteration algorithm to determine a 


more accurate solution. 


What is the general condition, in terms of g'(x), for fixed point iteration to converge 
to a solution of x = g(x). Given that the root of « = cos is approximately 0.739, do 
you expect a fixed point iteration algorithm to converge to the solution of x = cos x, 


or not? 


Re-write the equation x = g(x) in the form f(x) = 0 (explicitly including the cosx 
term). Then, write down the Newton-Raphson iteration formula for solving this specific 


problem. 


The following is an incomplete MATLAB script to determine the solution to the specific 
equation « = cosx using Newton-Raphson iteration. Complete the indicated lines 
of code so that the solution is determined, correct to 7 decimal places, and with a 


maximum of 20 iterations. 


xold = 0.75; 
for ii = ... %Incomplete line 1 
xnew = ... %Incomplete line 2 
if ... < 5e-8 %Incomplete line 3 
break 
end 
. = xnew; %Incomplete line 4 


end 
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Question 7 [10 Marks] 
Consider the following data: 
a OK ll 2 5 4 5 


te 0. 1, 26. - 3.6 5 7 
yi O 0.84 06 —0.44 —0.96 0.66 


(a) Draw a sketch that illustrates the approximate positions of the given data points in 


the xy-plane. 


(b) What is the minimum degree you would expect for a polynomial that passes exactly 


through all six data points? 


(c) If you were to estimate y(4) using linear interpolation, which pair of points would be 


the most appropriate to use? 


(d) Write down the Lagrange polynomial for x, = 5, that is, the polynomial L4(x) of 
smallest degree such that D4(a4) = 1 and Ly(2;) = 0 for all 7 € {0,1, 2,3, 5}. 


(ce) Assuming that L;(x), i = 0,...,5, are the Lagrange polynomials for the given x-values, 
write down the interpolating polynomial for the points given in the table as a linear 


combination of the polynomials L;(x). 


(f) Suppose, we want to construct a cubic spline f(a) through the given data points, that 
is, f is defined by f(x) = fi(x) for x € [x;, 2:41], 7 = 0,1,...,4, where the functions f,; 


are cubic polynomials which we can write in the form 
fi(x) = a(x — 24)? + b;(2 — 2)? + (x — 24) + dj. 


For the specific data given in the table above, write down the conditions on the coef- 


ficients a;, b;,¢;,d; that come from the constraints 
(i) f(z) = ys (ii) f3(va) = fa(a) 


(g) The MATLAB command, yy = spline(x,y,xx) returns a vector yy of values at the 
positions given in vector xx of a cubic spline interpolating the points (x(i) ,y(i)) for 
i=0O, 1, ..., n. Write a short MATLAB script to plot the spline over the interval 
(0, 10] and evaluate the spline at x = 4. 
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Question 8 [6 marks] From the values 


a 0 71/8 m/4 = 3n/8 a/2 


sin(z) O 0.3827 0.7071 0.9239 1 


a /2 
| sin x dx 
0 


using the trapezoidal rule and Simpson’s rule, both with 3 and then 5 function evaluations. 


compute approximations for 


Determine the exact value of the integral and compare it to your approximations. 


Please remember: This examination paper MUST BE HANDED IN. Failure to do 


so may result in the cancellation of all marks for this examination. Writing your 
name and number on the front will help us confirm that your paper has been 


returned. 


